Base pair dynamic assisted charge transport in DNA 
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An Id model with time-dependent random hopping is proposed to describe charge transport in 
DNA. It admits to investigate both diffusion of electrons and their tunneling between different sites 
in DNA. The tunneling appears to be strongly temperature-dependent. Observations of a strong 
(exponential) as well as a weak distance dependence of the charge transfer in DNA can be explained 
in the framework of our model. 
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Electronic transport is a ground to a wide range of 
important biological processes in DNA. Besides, the phe- 
nomenon has a fundamental physical interest, since the 
transport properties of biomolecules are expected to dif- 
fer considerably from those of macroscopic conductors. 
And at last, very recently material scientists also turned 
their attention to charge migration in DNA for the de- 
velopment of DNA-based molecular technologies. 

Although first attempts to measure DNA conductiv- 
ity Q , and to give a theory of the phenomenon ||^ have 
been made almost 40 years ago, the question concerning 
charge transport through DNA remains unsettled, and 
there is an impressive quantity of unexplained or par- 
tially explained data. Different publications report fre- 
quently contradictory results. Two kinds of techniques 
for getting information on charge transport in DNA are 
used. First, direct or indirect electrical conductivity mea- 
surements on micrometer-long DNA ropes are performed 
[^-0 . Experimental results obtained in this technique are 
ambiguous. DNA conductivity a was reported as almost 
metallic of the order il^^cm"^ [|j (in a recent publi- 
cation Q the authors claim they observed even proximity 
induced superconductivity in DNA) or semiconducting 
with a ~ 0.1 n^^cni^^ Very recently experimen- 

tal techniques have progressed to the point where the 
conductivity on individual lOnm long double stranded 
molecules was measured and the result implies that 
DNA is a good insulator. Clearly this frustrating situ- 
ation with conductivity measurements means that there 
are many relevant factors which can influence the charge 
transport in DNA by different way and which are hardly 
controlled in real experiments. The second technique, 
related to fluorescence quenching measurements on DNA 
strands, doped with donor and acceptor molecules [p|-^3[, 
seems more reliable, and it is our main concern here. 
In this technique photo-excitation of a donor associated 
with the stack of base pairs in some fashion, allows trans- 
fer of an electron to the stack. The migrating electron is 
trapped finally at the acceptor site, and charge transfer is 



monitoring by the yield of a chemical reaction accompa- 
nying the trapping process. The transfer rate is usually 
assumed to be characterized by a simple exponential law 
exp{— (3x) , where x is the donor-acceptor separation. Fit- 
ting to this law gives values of (3 ranging from 0.1 to 
I.4I-1. 

It is a common wisdom that DNA can be treated as an 
one-dimensional linear chain of stacked base pairs. We 
believe, that in the ground state every base pair con- 
tains bound electrons only. Then the charge is carried 
through DNA by excited electrons (or holes) which can 
jump between the base pairs. Below we propose a simple 
model which, to our meaning, reflects basic features of 
the electron transport in DNA. The picture includes the 
following ingredients: 

(i) The excited electron states at the base pairs are sep- 
arated by energy spacing larger than temperature and 
therefore thermally are practically not excited; 

(ii) Thermal motions of the DNA base pairs are elastic 
vibrations with a characteristic frequency ujb; 

(iii) Efficient charge transfer between neighboring base 
pairs takes place for rare events; 

(iv) The Coulomb interactions between the electrons and 
holes can be neglected for the description of hopping 
transport. 

Let us explain the point (iii) in more detail. For the 
static equilibrium DNA helix charge hopping is expected 
to be negligibly small since there is no significant elec- 
tronic overlapping between adjacent base pairs. Never- 
theless, sometimes, due to thermal fluctuations, exclu- 
sively favourable for hopping configurations of the base 
pairs occur, when an efficient hopping is possible. If the 
separation between the pairs is larger than the ampli- 
tude of their thermal vibrations, then probability of such 
events (which can be called "contacts") is small. Dura- 
tion of the contact can be estimated as the characteristic 
oscillation time uj^^. 

Now we discuss a correspondence of the assumptions 
enlisted above and experimental data. As a guide line we 
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use not only data known from the DNA literature but as 
well the data obtained for a wide range of organic linear 
chain polymers of stacked planar molecules (for a review 
see [Q). The contacts are related to mutual displace- 
ments and orientations of adjacent base pairs. Probably, 
hopping matrix elements are mostly sensitive to the rel- 
ative rotations (twist fluctuations) of the base pairs (see, 
e.g., Q). The characteristic frequency of these fluctu- 
ations, Wb, is usually estimated as being in the region 
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A small probability of the contacts is 
confirmed by experiment showing that the characteris- 
tic electronic hopping time t is larger than Wj^^, in the 
experiments ||,|^,|l||8| LUbT = 10^ 10^. Our first as- 
sumption (i) requires AE > T, where AE is the spacing 
in the spectrum of electron excitations for a base pair. 
The magnitude of AE can be measured directly, for the 
experiments [^-p3| AE > BOOK. Therefore the inequal- 
ity is satisfied. Rough macroscopical estimations of the 
Coulomb interaction Uc, as well as ab initio molecular 
orbital calculations of Uc, give few meV |26|, and thus 
Coulomb energy appears to be smaller (though of the 
same order) than hujb- We believe that it is enough to 
justify neglecting Coulomb interaction. 

The above reasoning leads to an Id hopping Hamilto- 
nian for the electrons 
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Here and aj' are electronic annihilation and creation 
operators at the site (i.e. the base pair) with the number 
i, and are the hopping amplitudes, which arc time- 
dependent quantities. The equations for the Heisenberg 
operators Ui are 
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We assume that different possess independent statis- 
tics, since are related to independent thermal pair base 
fluctuations. The hopping matrix element can be de- 
composed into a constant part (^i), that describes the 
coherent charge carrier motion in a completely rigid lat- 
tice, and a fluctuating part. Since the probability to jump 
is appreciable during rare events, the coherent part of ^ 
can be neglected in comparison with its fluctuating part. 

Note, that theoretical models based on hopping Hamil- 
tonians similar to Eq. (Q) are widely used to describe 
charge transport in solid state physics (see, e.g., p7| , p8|| ). 
For most of problems in this case the description cor- 
responds to electron migration in a steady energy land- 
scape, including thermally activated jumps over barriers 
and quantum tunneling through the barriers. It is quite 
different from our case. 

We assume that DNA molecules can be treated as 
homogeneous ones. Though the molecules are con- 
structed from four different nitrous bases, experimental 
data P |2^, as well as numerical first principle calcula- 
tions 29 3^ show, that the sequence of base pairs is 



not a decisive factor which determines electron transport 
in DNA. Quantitatively this condition can be formulated 
as 5E < hub, where 6E is an energy spacing between the 
(lowest excited) electron energy levels at different base 
pairs. The values of 6E, known mainly from numeri- 
cal electronic structure calculations |2^,|2^,^, are of or- 
der of meV. Thus SE is smaller than huJb for ujb given 
above, that justifies the picture. Besides some experi- 
ments (see, e.g., Q) are performed for artificial homo- 
geneous DNA, where 6E = 0. One expects that, due 
to the hopping, electron diffusion occurs on large time 
scales. For ^, treated as a white noise, it was demon- 
strated in the papers [^,^. Though our case is essen- 
tially different, there is good reason to believe, that the 
same behavior should be observed on time scales larger 
than the hopping time t. 

Below, we examine the particular case related to the 
fluorescence measurements, reported in the papers [p[-p3| . 
The donors are photo-excited and effects, related to the 
excited electron motion to the acceptors, are monitored. 
The energetic gaps SEd and 5Ea between the donor 
and the acceptor and the base pairs between them, are 
crucial for the hopping rate. The values of SEd and 
dEa (known mainly from ab initio numerical calculations 
[|l|,|6|j29[|3|]) can be estimated as lO^meT^. We see, that 
the inequalities SEd,dEa ^ huJb are satisfied. The elec- 
tron is always bound to the acceptor site more strongly 
than to a standard base pair, that is 6Ea > 0. As to the 
donors, the sign of SEd can be either positive or negative. 
If SEd is negative then a scheme of the electronic charge 
transfer from the donor to the acceptor is quite simple. 
Initially, the electron leaves the donor, jumping to the 
neighbor site, and then jumps between the standard base 
pairs, trapping finally at the acceptor. The case SEd > 
is more complicated. In order to have a driving force for 
the donor-acceptor charge transfer process the final state 
with the charge bound to the acceptor should be energet- 
ically favourable, that is the inequality SEa > SEd has 
to be satisfied. However, there are some base pairs in- 
between which play the role of the potential barrier for 
the electron. Therefore there are two possibilities for the 
electron to come to the acceptor. The first possibility is 
to jump initially from the donor to the neighbor site and 
then to move to the acceptor due to multistep hopping 
over the standard base pairs. The second possibility is 
the unistep (direct) quantum tunneling from the donor 
to the acceptor through the barrier. 

Since SEd 3> huJb, the probability for the electron to 
jump from the donor to the neighboring base pair due to 
dynamics of is negligible. At SEd > such a jump is 
possible if the electron absorbs a high-frequency phonon 
with the frequency LUph ~ ujd (= SEd/h). Correspond- 
ingly, at SEd < the electron jump from the donor is 
accompanied by emitting high frequency phonons. Such 
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dynamical vibrations with periods as short as tens fem- 
toseconds (i.e. phonons with ojph ^ lO^^s"^) were re- 
ported in the hterature ||3^-|35|]. Since htOph > T, oc- 
cupation numbers of such phonons are small. Thus, for 
6Ed > the probability for the electron to jump from the 
donor to the neighboring site contains two small factors: 
the probability of the contact and the probability to ab- 
sorb the high frequency phonon. It corresponds to the 
experimental situation where only a small fraction of the 
electrons are transported from the donor to the acceptor. 

When the electron leaves the donor, it starts to jump 
between the donor and the acceptor. It can return to the 
donor or can come to the acceptor. If 5 Ed < then the 
probability to return to the donor is negligible. We as- 
sume that even at 6Ed > the probability of the electron 
to jump to the donor or to the acceptor is smaller than 
the probability to jump to the standard base pair. There 
are two reasons for the assumption. First, the donors and 
the acceptors have chemical structures, different from the 
standard base pairs, that hinders for the contacts. Sec- 
ond, the jump has to be accompanied by the phonon 
emission, that diminishes its probability. The same is 
valid for the acceptor. Thus before being finally trapped 
at the acceptor site, the electron jumps many times back 
and forth over the base pairs between the donor and the 
acceptor, "smearing out" homogeneously over the all in- 
termediate base pairs. Then the relative probability for 
the electron to come to the acceptor is determined by 
the ratio of the probabilities for the electron to jump to 
the donor and to the acceptor from adjacent base pairs. 
This relative probability appears to be independent of 
the separation x between the donor and the acceptor. 
That explains why the rate of charge transfer sometimes 
is almost insensitive to the relative loading of donors and 
acceptors (see, e.g., ||3^ , |37|| ). The above picture imphes, 
that the total donor-acceptor charge transfer time should 
be larger than the electronic hopping time r, and it con- 
forms to experimental data (see, e.g., 

Now we consider the quantum tunneling for the elec- 
tron, strongly attached to the donor, that is the case 
SEd > 0. Though the potential barrier depends on time, 
at the condition SEd ^ hiVb the probability for the elec- 
tron tunneling from the donor to the acceptor can be 
calculated in the adiabatic approximation. To examine 
the tunneling process, one should consider the quasista- 
tionary electron state bound at the donor. In the spirit 
of our picture we assume ^ ^ cj^- Then the energy of 
the bound state is close to —SEd- Substituting dt by icod 
in Eq. (^, one obtains for the state 

Ui — ^-^ai^i if i > , a-i — — —ai^i if i < , (3) 

^d 

where we used the condition ^ <C cj^. Then the probabil- 
ity for the electron to be at the site n is determined by 
the average 
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Quantum averaging and averaging over statistics of ^ are 
performed at deriving Eq. (Q). Besides at the derivation 
we substituted (aq ao) « 1 justified by (a+a„) ^ 1. Note 
that the probability (Q) is determined by the simultane- 
ous statistics of ^. Remind that different are assumed 
to be statistically independent. Therefore the average in 
the right-hand side of Eq. (Q) is a product of 
For the standard pairs the quantities can be regarded as 
site independent ones. Therefore the probability of the 
electron to be at the nearest to the acceptor site is pro- 
portional to where n is the number of the 
standard pairs between the donor and the acceptor. 

The jump of the electron from the bound state to the 
acceptor is accompanied by the phonon emission. How- 
ever, the only x-dependent factor in the probability of 
the process is related to the average charge occupation 
number of a site n near the acceptor, established above. 
Thus, we obtain for the probability the exponential law 
exp(— /3a;) with 



/3 = a-iln(a-^/(|ep)) 



(5) 



Here ^ is the hopping probability for the standard base 
pairs and a = 3.4A is the distance between the base 
pairs in DNA. Note, that LOd depends on the donor type, 
whereas the average is mainly related to base pair 

vibrations. It follows from the above consideration, that 
the exponential law implies the condition (3 > . This 
conclusion is in agreement with the majority of pub- 
lished experimental data. Reported in the value 
j3 — 0.lA~^ (thus smaller than a~^) is, probably, related 
to an attempt to fit a complex behavior (including two 
processes: diffusion and tunneling) by a simple exponen- 
tial law. 

Let us stress that the quantum tunneling analyzed 
above is not a standard (static) tunneling described in 
textbooks. We have deal with dynamic tunneling which 
can be effective only when due to fluctuations of ^ there 
occurs some kind of a "bridge" from the donor to the 
acceptor. The exponential law, we found, is explained in 
fact by a small probability to have such a bridge, which 
is realized, when simultaneously many contacts between 
base pairs occur. In addition, the probability of this kind 
of tunneling is strongly dependent on the temperature 
via (ICP)- It is natural to assume that ^ exponentially 
depends on the relative displacement u of the neighbor- 
ing base pairs. Then (in the harmonic approximation) 
In(I^P) contains the term, proportional to (it^), which is 
proportional to the temperature T. Thus we arrive at 
the expression 



/3a = ci — C2T , 



(6) 
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where ci and C2 are temperature-independent factors. 
They can be extracted from the paper ci « 4, 

C2 « O.OIK^^. The values are in agreement with rough 
estimates ci ^ a/b, C2 ~ ks/ {Muj'^b'^), where ks is the 
Boltzmann constant, M is the base pair mass, and b is 
an electronic penetration length. It can be estimated as 
b ~ h/VmE ^ lA, where m is is the electron effective 
mass and E is its binding energy at the base pair. 

To conclude, for the electron, strongly bound to the 
donor, we established two different charge transfer mech- 
anisms: diffusion and tunneling. The diffusion leads to 
the charge transfer probability independent of the donor- 
acceptor distance x. However, the probability contains 
the small factor related to the electron jump from the 
donor to a neighboring site. The tunneling leads to the 
exponential dependence of the probability on x (with 
the temperature-dependent length /3"^). Therefore, it 
is not efficient for large distances. Thus, the exponential 
law has to be observed for small distances x whereas for 
large distances the charge transfer rate has to be inde- 
pendent of the donor-acceptor distance. Just this kind of 
behaviour was reported very recently ||2^. For the case 
of the electron weakly bound to the donor the hopping 
should always dominate over quantum tunneling. That 
explains why the rate of charge transfer sometimes does 
not behave exponentially even for small x |3^,0 . 

Note that in some cases the interaction of light, ion- 
izing radiation or chemically active reagents with DNA 
can result in loss of an electron at a specific site with 
formation of a hole. In this case the charge transport 
through DNA can be provided by holes (see, e.g., Q). 
The key issues for positive charge carrier transport are 
the same as for the electrons. As far as the physical pic- 
ture of charge transport is essentially the same for both 
kind of carriers, it can be described in the framework of 
the same approach. 
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